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Abstract 

. The problem of maximizing the average rate in a muhicast network subject to a coverage constraint (minimum 

. quality of service) is studied. Assuming the channel state information is available only at the receiver side and single 

CSl ' antenna nodes, the highest expected rate achievable by a random user in the network, called expected typical rate, 
is derived in two scenarios: hard coverage constraint and soft coverage constraint. In the first case, the coverage 

^! \ is expressed in terms of the outage probability, while in the second case, the expected rate should satisfy certain 

QQ . minimum requirement. It is shown that the optimum solution in both cases (achieving the highest expected typical 

■ rate for given coverage requirements) is achieved by an infinite layer superposition code for which the optimum 

^ ^ ' ' power allocation among the different layers is derived. For the MISO case, a suboptimal coding scheme is proposed, 

'"1 ' which is shown to be asymptotically optimal, when the number of transmit antennas grows at least logarithmically 

O . with the number of users in the network. 

'> \ I. Introduction 

00 ; 

^ \ The widespread application of wireless networks has motivated efficient transmission strategies for 



different applications. One of these applications is data multicasting where a group of users are interested 
^ in receiving the same signal, possibly at different levels of resolution. In a wireless multicast network, 



O 



X 



a common source is transmitted to users through a fading channel. In such networks, two criteria 
are usually studied as measures of performance: network coverage (minimum quality of service) and 
expected rate (typical quality of service). In the first criterion, the objective is to provide all the users 
with a minimum service regardless of their channel qualities. In the second criterion, the average data 
rate observed by a randomly selected user is considered where users with better channel conditions may 
receive higher data rates. An example for such networks is scalable video broadcasting in which all the 
subscribers expect to receive a basic signal, while those with better channel conditions might enjoy a 
higher resolution. 

In [1], the challenges in lossy multicasting are studied from an information theoretical point of view. 
In this work, for an analog Gaussian source with a bandwidth equal to the channel bandwidth, uncoded 
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transmission is shown to achieve the minimum average end-to-end distortion. The scenario in which 
the source has a larger bandwidth is studied in [2], where different methods of digital transmission are 
investigated. In [3], a different approach to broadcasting, called static broadcasting, is proposed. It is 
assumed that all the users receive the same amount of data from a common source, but with different 
number of channel uses as determined by their respective channel qualities. The actual transmission time 
in this scheme depends on the user with the worst channel, and hence, the transmission rate might be 
very low when the number of users is large. 

In this work, we consider a wireless multicast network in a quasi-static fading environment with additive 
Gaussian noise. The objective is to maximize the average performance, while a coverage constraint is 
satisfied. Average performance is defined as the Quality of Service (QoS) observed by a randomly chosen 
user (typical user), while the coverage requirement relates to the QoS observed by the user(s) with 
the worst channel condition(s). We assume that the transmission block is large enough to yield reliable 
communication. However, averaging over time is not possible because of the delay constraints. In other 
words, all the symbols within a transmission block experience the same channel gain. The channel state 
information (CSI) of each user is assumed to be known only at the corresponding receiver. For a channel 
with the above characteristics, the ergodic capacity is not defined, but the outage capacity is defined as the 
maximum rate decodable with a given probability [4]. In [5], a broadcast approach for a single user channel 
with these assumptions is proposed which optimizes the expected decodable rate. We apply "multicast 
outage capacity" and "expected minimum rate" definitions to characterize coverage in the network. Outage 
capacity is exploited when we have a hard coverage constraint. In this case, we require that with a given 
probability, within each transmission block, a specific amount of data is received by all the users. In the 
soft coverage constraint scenario, we relax the coverage constraint by stating it in terms of the expected 
minimum rate received by all the users within each block. For both hard and soft coverage constraints, 
another simultaneous criterion is the maximization of the expected typical rate which is defined as the 
average rate received by a randomly selected user. In general, there is a tradeoff between these two criteria. 

The minimum-service criterion has been studied in [6] for a single user fading channel, assuming CSI 
is known at the transmitter. In that work, given a service outage constraint for a real-time application, 
the average rate is maximized for a non real-time application sent on top of it. An adaptive variable 
rate code is proposed and shown to be optimum in that scenario. Similarly, a minimum rate constrained 
capacity measure is defined for broadcast channels in [7]. It is shown that the minimum rate region is 
the ergodic capacity region of a broadcast channel, with an effective noise determined by the minimum 
rate requirements. Using similar assumptions for the CSI availability, a more general case is considered 
in [8], where each user specifies its rate constraints in a triplet of maximum rate, minimum rate, and a 
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so-called shortage probability. 

In this paper, we use a broadcast model for an unknown quasi-static fading channel [5] and show the 
optimality of this model in our scenario. The same model is used in [9] to propose a multilevel approach 
for joint source-channel coding in a SISO channel (assuming the CSI is not available at the transmitter). 
The broadcast approach in [5] is extended in [10]-[12] for the case that the transmitter has partial channel 
state information, and in [13] for the case of MIMO block fading channel. [14] combines this broadcast 
approach with Hybrid Automatic Retransmission Request (HARQ) and shows that this approach results in 
high throughput and low latency in a point-to-point link. References [15], [16] use the broadcast approach 
in [5] in a two-hop relay network. 

We investigate the performance of the proposed scheme in both SISO and MISO cases. The MISO 
multicast asymptotic capacity limits are examined in [17], when the CSI is available at the transmitter. It 
is shown that the adverse effect of having a large number of users can be compensated by increasing the 
number of transmit antennas. In this work, we study a similar scenario and derive asymptotic capacity 
results for a large number of transmit antennas. 

The rest of this paper is organized as follows: In section 11, the system model is introduced. Section 
III focuses on the virtual broadcast model for an unknown fading channel when the network is delay 
limited. Sections IV and V discuss multicast networks when we have a single antenna at the transmitter 
and at each receiver. In section IV, we evaluate the optimum performance of the network in terms of the 
achievable tradeoff between the expected typical rate and the multicast outage capacity (hard coverage). 
Section V studies a similar problem of computing the optimum tradeoff, but for a soft coverage constraint 
where the expected minimum rate is used as the coverage criterion. Section VI investigates the MISO 
case, where we derive asymptotic capacity results for a large number of transmit antennas. Finally, section 
VII concludes the paper. 

Throughout this paper, we represent the norm of the vectors by ||.||, the conjugate transpose operation by 
{.y, and the expectation operation by E[.]. The notation "log" is used for the natural logarithm, and rates 
are expressed in nats. We denote fy{.) and Fy{.) as the probability density function and the cumulative 
density function of random variable y, respectively. Notation 1^ is used to define a binary function of x 
which is equal to 1 if event A occurs and otherwise. For given functions f{n) and g{n), f{n) = o{g{n)) 



is equivalent to lim„_ 



fin) 



0, and f{n) = uj{g{n)) is equivalent to lim„^oo = oo. We use A ^ B 
to denote the approximate equality between A and B, such that by substituting A by B the validity of 
the equations is not compromised. 
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II. System Model 

In this paper, we consider a wireless network broadcasting a common message. In the first part, it is 
assumed that a single-antenna transmitter transmits a common message to N single-antenna receivers. 
The received signal at the ith receiver, denoted by yi, can be written as 

yi = SiX + ni, (1) 

where x is the transmitted signal satisfying an average power constraint of E[x'^] < rii ^ CJ\f{0, 1) 
is the Additive White Gaussian Noise (AWGN), and Si ~ CN'{0, 1) is the channel coefficient from the 
transmitter to the ith receiver. The channel gain hi = which is assumed to be constant during a 
transmission block, has the following Cumulative Distribution Function (CDF): 

Fi{h)^F{h) = l-e-'*,Vi. 

The typical channel of the multicast network is defined as the channel of a randomly selected (typical) 
user. Since all the channels are i.i.d., the typical channel gain distribution satisfies 

Ftyp{h) = F{h) = 1 - e-\ (2) 

Since all the N channels are Gaussian and they receive a common signal, the coverage requirement is 
determined by the channel with the lowest gain hmui — minj(/ii), which is called the multicast channel. 
Due to the statistical independence of the channels, the gain of the multicast channel has the following 
distribution 

Pr |min(/ii) > /i} = {Pr {hi > h})^ = e-^^ 

As a result, we have 

Fmuiih) = 1 - Pr {min(/i,) > /i} = 1 - e-^^ 

In this paper, we deal with three quality measures defined as follows: 

• Multicast outage capacity, R^, is the rate decodable at the multicast channel with probability (1 — e). 

• Expected multicast rate, Rmui, is the average rate decodable at the multicast channel, i.e. Rmui = 

E[R{h)\h = h„iui], where R{h) is the decodable rate at the channel state h. 

• Expected typical rate. Rave-, is the average rate decodable by a randomly selected user, i.e. Rave — 

mm 
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III. Broadcast Model for an Unknown Quasi-Static Fading Channel 

In [5], it is shown that the expected rate for a receiver with a quasi-static block fading channel, unknown 
at the transmitter, and a stringent delay constraint, is equivalent to a weighted sum rate of a degraded 
broadcast channel with infinite number of virtual receivers, each corresponding to a realization of the 
channel. In this paper, we exploit the same model in a more general fashion. Noting our frequent use of 
this model, we first study it in more details. 

In this work, we assume a block fading channel for all users, where the channel state takes values 
according to a given probability density function (pdf) (assumed to be exponential) at the start of each 
block, stays unchanged during the coherence time (block length of the channel), and then changes 
independently at the start of the subsequent block. The channel state information for each channel is 
assumed to be available only at the corresponding receiver's side. In this case, if there were no constraint 
on the decoding delay, coding across different fading blocks would be possible, achieving the so-called 
ergodic capacity. However, in our model, we impose a decoding delay constraint which restricts the 
receiver to decode within a period equal to the length of a fading block. Each receiver decodes a fraction 
of the transmitted data which is supported by its corresponding channel. Hence, for any coding scheme, 
we have a function R{h) which determines the data rate decoded in channel state h. The time average of 
the rate decoded by a given receiver over infinite number of transmission blocks would be Efi[R{h)] [5]. 

Consider an infinite number of virtual receivers, such that receiver EXh is experiencing a fading level 
between h and h + dh. With these settings, EXh is receiving all the data received by EXh-dh, in addition 
to dRh, where 

dRh^ R(h) - R(h-dh). 

The virtual receivers introduce a degraded virtual broadcast network in which the rate associated with user 
RXh is dRh. From the degraded nature of the Gaussian broadcast channel [18], it follows that dRh > 0. 
The original receiver corresponds to the virtual receiver RX^ with probability r){h)dh, where ri{h) is the 
channel gain probability distribution function. 

With this interpretation, for a given coding scheme, the distinction between different channels introduced 
in the previous section, namely multicast channel and typical channel, is translated to the difference 
between the probability distribution based on which the original receiver is represented by the virtual 
receivers. Note that both the multicast and typical channels correspond to the same common virtual 



broadcast network, and the measures defined in the previous section could be written as follows B 

/■oo 

Rave = / Rih)7]{h)dh 

Jo 

poo 

= -il-Fih))Rih)\'^+ (1 - F(/i))di?, 

(1 - F{h))dRH, (3) 



(a) 







oo 



R^ui = / {1 - Fmui{h))dRh, (4) 



R, = R{K)= [ ^ dRh, (5) 
Jo 

where = F~li{e) = — . in the above equation, (a) follows from the facts that R{0) = and 

lim/,^oo(l — F{h))R{h) = 0, where the latter is due to the fact that 1 — F(h) = and R{h) grows 
logarithmically with h. A similar argument is applied to conclude @i. In the case of the multicast channel 
in ([5]), the original receiver has a channel gain less than with probability e, and hence, the highest 
decodable rate is R{h^), with probability 1 — e. 

As seen above, the performance measures in this work are different weighted sum rates of the virtual 
broadcast network, resulting in vectors {R^, Rave) and (Rmui, Rave) for the hard and soft coverage con- 
straint scenarios, respectively. We refer to the maximum achievable regions for these vectors as the capacity 
region corresponding to each case. Next, we provide a definition for the optimality of a performance vector. 

Definition 1 The boundary set B of a closed convex region M C M" is defined as 

= {x G e X + X G ^} (6) 

where and M_|_+ are the set of nonnegative and strictly positive real numbers, respectively. 

With the above definition, a performance vector is optimal if it is in the boundary set of all possible 
performance vectors. In the following theorem, we show that the optimal performance vector for each 
coverage constraint scenario is achieved by superposition coding in which the rate of the virtual receiver 
RXh is given by 

d/?. = log|'l + -^^#^l, (7) 



1 + 



u 



where p(.) is the power distribution function such that p(h)dh is the amount of power allocated to the 
virtual receiver corresponding to the channel state h. 



'with a small misuse of notation, to simplify the problem formulation, the integration in the right hand side of all of the three equations 
is computed over h, where dRh is expressed as an explicit function of h throughout the paper. 
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Fig. 1. A schematic figure for Lemma 1. 

Theorem 1 The boundary set of {W\, W4), where Wi = Wi{h)dRh Vz /i' a positive weighted sum 
rate of the underlying virtual broadcast channel, is achievable by a super-position coding scheme. 

Proof: To prove the theorem, we first state and prove the following lemma: 

Lemma 1 Consider a mapping function g{.) from a closed region ^\ C to C M*^ such that 
glx.) = Mx, where M G x M", such that each row of M contains at least one positive element. 
Denote Bi and B2 as the boundary sets of the regions and ^2, respectively. Then, we have 

B2 C g{Br) 

Proof: Assume this is not true. Hence, there must exist X2 G B2 such that X2 ^ g{Bi) and xi G ^i, 
such that X2 = (/(xi). Since xi ^ Si, there exists x'^^ G Bi such that x'^ — xi G Defining Xg = 

g{^i) G =^2, we have 

x'2 - X2 = M(x; - xi) C M^,,, (8) 

which contradicts the fact that X2 is in the boundary set of ^2 and the lemma is proved. ■ 
Figure ^ illustrates an example of Lemma [H when n = 2 and k = \. The region ^1 in this example 
is defined by xi,X2 >Q,x\ + x\ < 1 and the mapping is defined by g{xi,X2) = Xi + X2. The boundary 
of ^i, i.e., Bi, is the quarter-circle xi,X2 > 0, xl + X2 = 1, which is mapped to the line segment [1, ^/2]. 
The boundary of M2 = g{^i), i.e., B2, which is the line segment [0,-\/2], is the point \/2, which is a 
subset of g{Bi). 

We can arrange the set of the rates of the virtual broadcast network in a rate vector dR = [d_R/i]jf and 
use the result of Lemma [H when n tends to infinity. In this case, the matrix transform will tend to k 
weighted sums of the infinite dimensional vector x = dR as follows: 



Mx 



Wi{h)dRh 



Wk{h)dRh 
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Setting x{h) = dR^, we can conclude that the boundary region of (Wi,- ■■ ,Wk) is a subset of the 
transformation of the boundary set of the rate vector dR, which is achieved, as shown in [19], by 
superposition coding. In other words, any vector v in the boundary set of {Wi, ■ ■ ■ , Wk), is achieved 
by a multi-layer code in which the rate of the virtual receiver HXh is equal to 

hpv{h)dh \ 
^ + hQ^^p^{u)duJ ' 

where the power distribution function pv{h), indexed with v to emphasize its dependence on this vector, 
satisfies py,{u)du = . This completes the proof of Theorem [TJ ■ 
Using the above theorem, it easily follows that the optimal performance vectors [R^, Rave) and {Rmuh Rave), 
defined in the sections IV and V, respectively, can be achieved using superposition coding. 

IV. Hard Coverage Constraint 

In this section, we consider a scenario where a minimum rate of R^ (multicast outage capacity) should 
be delivered to all the users with a probability of (1 — e), where e is the probability of outage. Given this 
constraint, we want to maximize the average rate received by a randomly chosen user in the network, i.e. 
Rave- This includes the rate received by a typical user, even if the user is in outage. We can categorize 
the possible states of the virtual broadcast network in two group: (i) h < h^, in which case the rate dRh 
contributes to both R^ and Rave, and (ii) h > h^, in which case the rate dRh contributes only to Rave- 

The problem of maximizing Rave without any constraint on R^ is studied in [5]. On the other hand, 
from [4], we know that the maximum R^ without any constraint on Rave is achieved by a single layer 
code with power ^ and rate 

a = l0g(l + /ie^). (9) 

In this section, we investigate the tradeoff between R^ and Rave- 
Setting wi{h) = l{h<h^} and W2{h) = 1 — F{h), Theorem [U states that the boundary set of {R^, Rave) 
is achieved by superposition coding, in which 

dR^, =los(l + = i-m+pih)dh Mp_ 

and I{h) = /^^^ p{u)du. Note that dRh is not necessarily very small since our power allocation function 
might have some impulses. As stated earlier, we want to jointly optimize the weighted sum of these rates 
according to the weighting functions wi{h) and W2{h). The optimization is on the function p{h), however, 
in the following we show that it can be simplified to a point optimization problem. 
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Definition 2 The channel gain-interference function, s{p) is defined as 

s{p) = sup{h\ I{h)>p}. (11) 

In fact, the channel gain-interference function s{p) is the inverse of the interference function in terms 
of the channel gain in the points it is invertible and determines the channel gain of the virtual receiver 
experiencing the interference level p. It is evident that s{p) is a decreasing function of p. According to 
(flOl ). we can write the expected typical rate as 



Ra., = - F{h))dRh = f^^ g{p, s{p))dp, (12) 

where g{x,y) = (1 — F{y))j^^. In deriving (fT2l) . we have used the fact that the contribution of the 
virtual receiver RXh in the expected typical rate, from ([3]), can be written as (1 — F{h))dRh = (1 — 
^(^)) filh)^''^^^'^^ Trk- Noting that in the interval p E [I{h), I{h)+p{h)dh], h can be written as s(p)B, the 
contribution of RX^ in Rave can be written as jllj^-^^''^^^^'^ 9{p, s{p))dp. Since Rave is the summation of the 
contributions of all virtual receivers, we have Rave = Ii(h)^^^^^'^'^ diPi s{p))dp = //(g^"* g{p, s{p))dp = 
if 9{Pi s{p))dp, which is due to the fact that /(O) = ^ and I{oo) = 0. 
Differentiating g{x,y) with respect to y, we obtain 

In the case of Rayleigh fading, we have f(y) = 1 — F(y) = e^^ . By studying the behavior of -§^g{x, y), 
it is easy to show that -i-g{x,y)\x=p is positive for y < Iq^{p) and is negative for y > Io^{p), where 



Uh) = ^^-T^r.^^'^ - As a resuk. 



fjh) , 

h^f(h) 

arg max{g{x, y)\x=p) = lo^ip)- (14) 

Note that /o(.) S is indeed the interference function corresponding to the optimal power allocation in the 
unconstrained maximization of Rave solved in [5]. 

Definition 3 The multicast interference level a for a given channel gain-interference function s{.) is 
defined as 

a = min {p\s{p) < /i^} • (15) 

In fact, a is the level of the interference observed by all the virtual receivers contributing to R^ due to 
the power allocated to the upper levels not contributing to R^. It follows from definitions [2l and [3] that 

^Note that this is true even if p(.) contains an impulse at h. 

^Note that the function /o(0 is monotonically decreasing and as a result, it is invertible. 
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a = Using this definition and the same arguments as in (fT2l) . in (|5]) can be written as 



R. = I dRh 







h, i-i(h)+p(h)dh f^^p 
Jiih) 1 + hp 



Ji{h) 1 + s{p)p 

^(''^^ s{p)dp 
/(o) 1 + s{p)p 

m{p, s{p))dp, (16) 



a 







l{s(p)<h,}m(p, s(p))dp (17) 



where 



y 

'm{x,y) = — , (18) 

1 + xy 



and (flTl) follows from the definition of a in (|T5I) . Here, we assume that < 1. This is not a restricting 
assumption as the solution in the other case {h^ > 1) follows using a similar approach. However, this 
assumption simplifies the derivations as it guarantees loiK) > 0. On the other hand, this assumption 
is equivalent to e < 1 — e^^, which covers most of the cases of interest since we expect the outage 
probability to be much smaller than 1. Using (fT2)) and (fT6l) . the problem is translated to 

max s{.) Rave = niax^(.) jf g{p, s{p))dp, (19) 

subject to 

R,= / m{p,s{p))dp>Ca, (20) 



where is defined in and is a normalization factor which expresses R^ in terms of C^. Since it 
is not possible to achieve values of R^ above C^, we can restrict ourselves to the values of ( between 
and 1. Note that the maximization in (fT9l) is over all decreasing positive functions s(.). Also note that a 
in (|20l) is the multicast interference level defined in (fT5l) and depends on the function s(.). 

From (fTSi) . we note that for any chosen x, m{x, y) is an increasing function of y. Hence, noting the 
definition of a in (fT5l). we can write 



R^< mip,K)dp = log^T^J 



-log(l + /i,a). (21) 
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Therefore, following (|20l ) and (|2TI) . we have 

CCe < a - l0g(l + /lea) 

p(i-C)a _ 1 

< • (22) 

he 

Lemma 2 For the optimizer of rfiPl). we /zave a < Io{K). 

Proof: Assume a > Io{h^). Denote the optimizer function by s*(.) and its resulting multicast outage 
capacity and expected typical rate by R* and -R*„e, respectively. Also, define s{p) as 

' l,\p) p<h{K) 

sip) = I K Io{K) <P<a , (23) 
s*{p) p > a 

and its resulting multicast outage capacity and expected typical rate by and Rave, respectively. We can 
write 

rCe — JrC^ — 

[l{s(p)<h,}m{p,s{p)) - l{s'{p)<h,}m{p,s*{p))] dp = 
iTi{p, s{p))dp > 0, (24) 







a 



where (a) follows from the facts that (i) s{j)) < for p > Io{h^), (ii) s*{p) < for p > a (from 
the definition of a in (fT5])). and (iii) for p > a, s{p) = s*{p) (from (|23])). (|24l) implies that R^ > R*^. 
Moreover, we have 

Rave-Kve = / id {P , Kp)) ^ 9 (P , S* (p))) dp (25) 
^0 



(a) 



[dip, lo^ip)) - 9iP, s*{p))] dp 



+ / {9{P,K) - 9{p,s*{p)))dp > 0, 

where (a) follows from the fact that for p> a, s{p) = s*{p). In the above inequality, the positivity of the 
first term in the left hand side is concluded from (fT4)) . and the positivity of the second term is concluded 
from the fact that s*(p) > for p < a (from the definition of a in (fTSi) ). and also (fT3l) which implies 
that g{x,y)\x=p is decreasing for y > Iq^{p). Therefore, Rave > Rive ^^'^ -^e > which contradict our 
assumption of optimality of s*(.) and the lemma is proved. ■ 
The above lemma states the fact that, applying the multicast outage constraint, more power will be 
allocated to the channel gains lower than the outage threshold, compared to the unconstrained scenario 
studied in [5]. 
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Lemma 3 Given a, the optimizer of dTPI) is given by 



K a<P<ltihe) > (26) 

it~\p) P>Ix{K) 



where Ix^{h) = ^^^^ hy(h) ^'^''^^ which the superscript {•Y''^ stands for the "hard coverage" constraint 
scenario, I^'^ ^(.) represents the inverse of the function I^^i-), and 



A 



0, if m{p,Sa{p))dp> (Ce 

arg (J^ to(p, Sc<(p))dp = CCe), otherwise 

Proof: As observed from (|20l) . the value of s{p) in the range < p < a dose not affect the mukicast 
constraint. Hence, (fT9l) can be written as 

pa 

max Rave = max / g{p,s{p))dp (27) 

„^ 

+ max / s(p))dp 

s{.)jfmip,s{p))dp>(:cJa 

Note that as mentioned earlier, all the maximizations are performed over positive decreasing functions. 
Moreover, since s{a) = (due to the definition of a in (fT5l)). the solution of the above maximization 
problem must satisfy s{p) > in the interval p E [0,a] and s{p) < elsewhere. (fT4l) implies that the 
first term in ([27] ) can be upper-bounded as follows: 

max / g{p,s{p))dp< / ^((p, /o"^(p))dp. (28) 
<■) Jo Jo 



Moreover, writing K.K.T. condition for the second term of (|27|) . the maximization problem will be 
translated to 

max / Tx{p,s{p))dp, (29) 

<■) Ja 



where Tx{x,y) = g{x,y) + Xm{x,y). A is 0, if the outage constraint is not limiting; otherwise, its value 
can be computed using the ou 
with respect to y, we obtain 



can be computed using the outage constraint m{p^ s(p))dp = C,C^. Differentiating the function Tx{x, y) 



d^. . X + l-Fiy)-yf{y)il + xy) 

In the case of Rayleigh fading, by studying the behavior of ^Ta(x, y\x=p), it follows that ■^Tx{x, y\x=p) > 
OfoTy< Ii^'\p) and ^Ta(x, yl=p) < for y > If "'(p), where (/i) = As a resuh, 

arg max TA(a;,j/)U=p = "^(p). (31) 

''Note that the function is monotonically decreasing for h < he (which is the region of interest here), and as a resuh, it is invertible. 
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Note that as Io{h) < I\{h), yh, and a < Io{h^) (from Lemma O, it follows that a < I^^(h^). Hence, 
the integral Tx{p, s{p))dp can be written as the summation of the two integrals jj^ ^'^^'^ Tx{p, s{p))dp 
and )Tx{p, s{p))dp. First, we note that in the whole interval of [a, we have s{p) < (due 

to the definition of a). Moreover, as /^^ ^(.) is a decreasing function (in the interval [a, ^]), we have 
he < ^(p), in the interval [a, /^"^(/ij]. Combining this with the fact that s{p) < in [a,I^^{he)\ 
implies that s{p) < I^'^ ^{p) in this interval. As seen before, Tx{x,y) is an increasing function of y for 
y < ^(x). Consequently, 



Tx{p,s{p))dp < / Tx{p,K)dp. (32) 



Moreover, using (|3TI) . we have 



/ Tx{p,s{p))dp< / Tx(p,I^'~\p))dp. (33) 



Hence, for any function s(p) such that s{p) < for p > a, we can write 



I Tx{p,s{p))dp < / Tx{p,Sa{p))dp. (34) 



Combining the above equation with (1281) yields that 



max < / g{p,Sc,{p))dp. (35) 
<■) Jo 

To complete the proof, it is sufficient to show that Sq,(.) satisfies the conditions mentioned earlier; i.e., 
Sa{.) is a decreasing function and Sa{(y) = h^. The latter is obvious from the definition of Sq,(.) in (|26l) . 
For showing the former, we first note that Iq^{-) and Ix'^ ^(.) are decreasing functions in the intervals 
[0,a] and [Ix'^{he), respectively. Moreover, from Lemma [2l we have a < /o(/ie) which implies that 
Io^i<^) > K- This shows that Sq(.) is decreasing in the whole interval [0, which completes the proof 
of the lemma. ■ 
An interesting consequence of Lemma [3] is that the problem of maximization over the function s(.) is 
simplified to the point optimization problem over the value of a. 

Theorem 2 The capacity region of a Rayleigh fading multicast network with a hard coverage constraint 
is given by 

Rave < max / g{p,Sa{p))dp, 

(36) 

Re < Ca, (37) 
where ( changes from to 1. 
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Proof: The proof follows from Lemma [2l Lemma [3l and inequality (|22l) . ■ 

Corollary 1 For any outage probability e > such that he < Iq^{^), the capacity region of a Rayleigh 
fading multicast network, i.e., (i?^, Rave), is given by 

Re < ^og(l + —^^4^—-], (38) 



Rave < 2{EMP)) - m)) - (e-'^^^ - e-') 

+ e-'^\og(l + —^^i^—-), (39) 



where B changes from to 1, 6(3) = -, — 2 and EUx) = ('°° ^dt. 

Proof: Since h, < Iq\^), it follows that Iq{K) > Noting that /o(/i) < /a'(/i), V/i, A > 0, it 
follows that Ix%K) > ^ for any A > 0. Therefore, (|26l ) can be written as 

. . J Io\p) P<o^ 

= < . (40) 

I /i, a<p< ^ 

In this case, the solution to the maximization problem (l36l) is a = 77— ^- Defining /? = 1 — 
first we note that < /3 < 1. Moreover, (|37l) can be written as 



Re < CCe 



a-iog(i + /i.(i- 

logfl+ , , )■ (41) 



l + /ie(l-/5) 

The interference function corresponding to Sa{p) in (|40|) can be expressed as 

/(/i) = aU {Iq^ (a) - /i) + 

-h) + Io{h)U [h - I^' (a)) , (42) 

where f/(.) denotes the unit step function and a = (1 — P)^. Differentiating I{h) with respect to h 
results in the following power allocation function: 

p{h) = {l^-a)6{h-K)+poih), 

where 



Po{h) 



^-^ I,\a)<h<l 
else 
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is the power allocation function in the unconstraint problem [5]. Using ([3]) and (flOl) . the expected typical 
rate can be written as 



/■oo 

Rave = / {l-F{h))dRh 

Jo 



-h , A \ 

io\a) 1 + hIo{h) 

2— d/i- / e-'^d/i. (43) 

Noting that I^^ia) = 6{a) = i^^^^:^ and a = (1-/5)^ completes the proof of the corollary. ■ 
An interesting conclusion of Corollary \\\ is that, the expected typical rate is maximized when the 
multicast minimum rate is provided in a single layer. In the case we have no multicast constraint, it is 
shown in [5] that a multi-layer code with a small rate in each layer is optimal in terms of maximizing 
the expected rate. However, when we are constrained to distribute a fraction of the available power to a 
set of low channel gains [0, h^] (coverage constraint), the optimum solution allocates all the power to the 
highest gain, i.e. h^. 

Note that the assumption < Iq^{^) is not difficult to satisfy, since the outage probability e is usually 
small. Moreover, as = — , the value of decreases significantly with the number of users. For 

example, for = 5 and ^ = 100, the outage probability e could be as high as 0.38 in order to have 

< Figure (O shows the capacity region of this network when e = 0.01. It is evident that due 

to the hard coverage constraint, the multicast outage capacity is in general very small in comparison with 
the expected typical rate. 

V. Soft Coverage Constraint 

In the previous section, we observed that a strict coverage constraint can result in very small values for 
the multicast outage capacity. We can relax the coverage requirement by relying on the expected multicast 
rate, i.e. Rmui- This results in the performance vector {Rmuh Rave) and its corresponding capacity region. 
According to Theorem [H the optimality of superposition coding is concluded for this performance vector. 
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0.05 0.1 0.15 

R (nals/symbol) 



Fig. 2. Hard coverage constraint: multicast outage capacity vs. expected typical rate for ^ = 100 and N — 5. 



Theorem 3 The capacity region of a Rayleigh fading multicast network with soft coverage constraint is 
given by 



Rn 



R 



mul 



1 + uP^^{u) 

-Nu Up'^^Wu 

1 + uPAu) ' 



(44) 
(45) 



where 



ifh<hQ 



e-fe(l-fe)+^e-iVh(l-Arfe) , , , 

h> hi 



(46) 
, and 



in which the superscript {.y^ stands for the "soft coverage" constraint scenario, pt^{h) — 
ho and hi are real numbers, such that 

e-^o(l-feo)+7e-^^o(l-Affeo) ^ ^ 
/i2(e~''0+7Are-'V''o) 

e-''i(l-fei)+7e-^'"i(l--Affei) _ p, 

for different values of •y > 0. 

Proof: If we set wi(h) = 1 — Ftyp{h) and W2{h) = 1 — Fmui{h), Theorem [T] states that we should 
search for an infinite layer superposition code. Considering p(h)dh as the power of the layer associated 
with the channel gain h, the corresponding rate ijf| 



1 fi_^J^PWdh\ hp{h)dh 



(47) 



Note that here, unlike the hard coverage constraint scenario, the power distribution function does not have any impulses. In fact, the 
optimization problem in the soft coverage constraint scenario, as seen in the proof, is similar to that of the unconstrained scenario [5] for 
which the optimal power allocation function has been proved to have no impulses. 
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where 

PCX) 

I{h) = I dup{u), 
Jh 

and 

/(O) = ^. 

Using the above expression, the rate corresponding to the fading level h is 

up{u)du 



R{h) 

Following the definitions, we have 



1 + ul(u) 



oo 



Nu up{u)du 



Rmui= {l-F^^i{u))dR{u)= e-^ \ ' (48) 



JO 

oo POO 



1 + ul{u) 



Ra.e= I {l~Fty,{u))dR{u)= I e-^-^^^^rr- (49) 
Jo Jo 1 + 

The problem is that given Rmui = r, what is the maximum achievable Rave- In other words, 

fl„. = m»re-"i5^, (50) 

/(«) Jq 1 + ul{u) 



subject to: 



^_,„u^ (51) 
1 + ul[u) 



/(O) = ^, 

and 

/(oo) = 0. 



Equivalent 
problem: |^ 



y, to derive the capacity region [Rmuu Rave), it is sufficient to solve the following optimization 



ra&-KRave + iRmuh (52) 
!{■) 



subject to 

/(O) = ^, 
/(oo) = 0, 



*Note that as we are allowed to user time-sharing, the capacity region is convex. As a result, the capacity region can be characterized as 
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Fig. 3. Soft coverage constraint: expected minimum rate vs. expected typical rate for ^ — 100 and N — 5 



for all values of 7 > 0. To solve the above optimization problem, we define S{x, I{x), I'{x), 7) as follows: 

S{x,I{x),I'{x),-f) = e~ \ ""^'^f. , +76-^" ^^'^^^ 



1 + xl{x) 1 + xl{x) 

(53) 

Note that 

I'{x) = —p{x). 

The necessary condition for I{x) to maximize (l50l) with the constraint (|5T1) is the zero functional variation 
[20] of S{x,I{x),r{x),^), 

d ^ d ^ _ ^ ^^^^ 

dl dx dl' 

where 

dx dl' ~ l+xl{x) ^ ) (l+x7(x))2 • 

Therefore, (|54|) simplifies to a linear equation which leads to the optimum interference function given in 

■ 

Figure ([3]) shows the capacity region for = 5 and ^ = 100. It is observed that the maximum Rave 
is achieved for Rmui < 105. It is shown in [21] that a good fraction of the highest expected rate with 
infinite layers is achieved by two layers. Figure ([3]) shows that this is true for our multicast network as 
well. Furthermore, we observe that the two-layer rate region gets closer to the capacity region for higher 
values of 
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Fig. 4. Power distribution function p{h) for no multicast requirement and for Rmui = 1-4 



Without the multicast constraint, the power distribution function to maximize the average rate would 
be [5], 

j w-^ ^o<h<l 
Po{h) = < 

I else 

where 

2 

So 



1 + VI + 4^ 

This function is depicted in figure (0]), and is compared with the case corresponding to a multicast 
requirement of Rmui = 1-4. As shown in the figure, the coverage requirement has shifted the power 
distribution to smaller channel gains to provide service to users with worse channel conditions. 

VI. Extension to MISO 

In the case that there are multiple (M > 1) antennas at the transmitter, one can adopt the broadcast 
approach proposed in [5]. In this approach, the receiver with unknown quasi-static fading MISO channel, 
denoted by h, is modeled as a continuum of receivers each associated with a channel realization. These 
receivers are ordered in a degraded fashion. However, since MIMO Broadcast Channel (MIMO-BC) is 
inherently non-degraded [22], this approach dose not necessarily lead to the optimum performance. 

Assuming a single antenna at each receiver, the ordering of the virtual receivers in this approach is 
based on their normalized channel norm, i.e., ^^^g^. Hence, the rate of the virtual receiver at the fading 
level equals 

ihht||\ , 



M I iug, ^ 1 -r ^ ^ p^, llhht|| 



where Ps and Pj are the decodable and undecodable signal power levels, respectively. 
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Now, assume there are N users in the network, all receiving a common source through an infinite-layer 
code. Same as before, we would like to design the code to maximize the average rate observed by a 
typical user, while providing a given coverage constraint for all the users. For this purpose, we should 
provide the minimum rate to the worst user in the degraded broadcast model, i.e., the user with the lowest 
channel norm. The normalized channel norm of user i, denoted by is a scaled random variable 

with 2M degrees of freedom, with the following CDF: 



Ftyp{h) 



llhht 



ih) 



r(M, Mh) 

r(M) ' 



(55) 



where r(.) is the Gamma function, and r(., .) is the upper incomplete Gamma function [23]. Since the 
users' channels are statistically independent, the distribution of the norm of the worst channel can be 
computed as 

p J ■ llh^hlll . 
rr < mm — > n 

i M 

rT{M,Mh)\ 

V r(M) ) 




Hence, the CDF for the worst channel norm is 

Fmul{h) — 1 



r(M, Mh) 

r(M) 



N 



(56) 



Here, we just consider the soft coverage constraint scenario. Following the same approach as in section 
rv, we obtain 



■mul 



[1 - Frr,ul{u))dR{u) 



r(M, Mu) 

r(M) 



up{u)du 
1 + ul(u) ' 



(57) 



Rn 



(1 - Ftyp{u))dR{u) 
T{M,Mu) up{u)du 



(58) 



,0 r(M) 1 + 

where p{u) and I{u) are the corresponding power allocation and interference power functions, respectively. 
To characterize the achievable rate region {Rave-, Rmui)', we need to solve the following optimization 
problem for all 7 > 0: 



maxima 

!{■) 



iR 



muh 



(59) 



subject to 



m 



/(oo) 
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Defining ^(x, /(x), /'(x), 7) as 



S'(x,/(x),/'(x),7) 



T{M,Mx) xl'ix) 



r(M) 1 + xl{x) 



+ 7 



r{M,Mx)\^ xl'ix) 



N 



T{M) J 1 + xl{x) 

and setting its functional variation equal to zero to maximize the average rate, and defining 

^ T{M,Mx) f r{M, Mx) \ 

r(M) +H r(M) J 



w [x) + 



similar to (1541) . we obtain 



which implies that 



1 + xl{x 



I(x) 



1 +XJ(X))2 

w{x) 1 



x'^w'{x) X 



Noting that -^T{M,x) = —x^ ^, we can write 



w'{x) 



r(M) 



1 + 7A^ 



r(M, Mx) 

r(M) 



Substituting in (|62l) yields the optimizer I{h) as 



^ if /i< /lo 

ho < h < hi 
/i > /ii 



where 



r(M, M/i) +7^^ 



r(M,Mfe) 



(M) 



JV 
JV-l 



and /lo, /ii are the solutions of the following equations: 



N-l 



1 

/I' 



/i(/io) 

^{hi) 



0, 



(60) 



(61) 



(62) 



(63) 



(64) 



respectively. 

The achievable rate region is shown in figure ([5]) for different values of N and for M = 2. From this 
figure, we can see that as the number of users decreases, the proposed achievable rate region expands. It 
is also evident by comparing the regions for the MISO and SISO cases with = 5 users (figures (|3]) and 
©) that using multiple antennas improves the achievable rate. However, its effect is more considerable for 
the multicast channel as compared to the typical channel. This prominent gain for the multicast channel 
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Fig. 5. Soft coverage constraint: MISO expected minimum rate vs. expected typical rate for different number of users, M = 2 and ^ — 100 

is due to the fact that we are using multiple independent paths to convey data (higher diversity order), 
so the probability of having very low channel gains significantly decreases. In fact, we will show that 
one can compensate the effect of having a large number of users by increasing the number of transmit 
antennas. More specifically, if both and M tend to infinity and M grows fast enough with respect to 
A^, we show in the next theorem that the multicast rate can reach the expected typical rate and our scheme 
gives the optimal solution. 

Theorem 4 For large values of M and N, the proposed infinite layer superposition coding will provide 
Rmui such that 



where Ropt is the highest achievable average rate for a typical user in the network, a is an arbitrarily 
small positive number, and denotes any function of N which tends to infinity as N oo. 

Proof: First, we derive an upper bound on the achievable average rate for a typical user by assuming 
no stringent delay constraint, meaning that the transmission block can be chosen as large as the fading 
block. In this case, the channel has an ergodic behavior, and hence, the ergodic capacity is defined and 
is shown to be 




(65) 



if 



M > 



(1 + ^)2^2 



(66) 




(67) 



As a result. 




(68) 
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Using the concavity of log function, and having the fact that E 

Ce., <log(l + ^). 



||hht| 
M 



1, we have 



(69) 



We will show that our scheme provides a multicast rate arbitrarily close to this upper bound, if we use 
enough transmit antennas. Since this upper bound is larger than the expected typical rate, the theorem 
will be proved. Using a single-layer code with power ^ [] and rate R„, where 



i?^ = log(l + ^(l-a')), 



(70) 



and 



a 



the expected multicast rate in our network will be 



R 



mul 



Pr 



M 



>l-a'} 



(71) 



where ||hh^||m^i; = minj ||hjh|||. Regarding the central limit theorem [23], the distribution of where 



TTllhh^l 



/ll^ + /l2^ + ■•• + 

M 



(72) 



and {/im}m=i's independent Rayleigh distributions with unit variance and unit mean, approaches to a 
Gaussian distribution with the CDF 

' h-1 



F{h) = Q 



Consequently, the CDF of the multicast channel will be 



(73) 



M 



h-l 



N 



M 



Fmul{h) ^ 1 - Q 

Using the above equation, dTTI ) can be written as 

R^ui ~ Q{-^a'fR„ 



(74) 



N 



Ra- 



(75) 



1 - QiVMa') 

Assuming M is large enough to have \jMci' >> 1, and consequently Q{-\/Ma') << 1, we can rewrite 
the above equation as 

^Note that as the single-layer coding is a special case of superposition coding, our proposed scheme outperforms the single-layer coding 
scheme. 
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Now, using the approximation 



for large values of x, we can write 



Q{x) ^ -^e-'r (77) 



Therefore, having 



Q{VMa')<e — —. (78) 



2\ogiN)+coil] 



(7'2 



M = ^, (79) 

incurs NQ{-\fMa') = o(l), and as a result. 



Urn Rr^^i -R„ = 0. (80) 

AT— >oo 



Moreover, assuming a <^ 1, (|70|) can be written as. 



R„ ^ log(l + <^) 



> Cerg-^, (81) 



where the second line results from (|69l) . Combining (I68I) . (|80|) . and (fSTI) . the result of Theorem 4 easily 
follows. ■ 
Theorem |4] simply implies that as the number of transmit antennas grows at least logarithmically with 
the number of users, the gain of the worst channel converges to the gain of the typical channel in the 
network, with probability one. In other words, increasing the number of transmit antennas provides fairness 
in the system, such that all users almost get the same quality of service. This fact is also noticed by [24] 
in the context of MIMO-BC. 

VIL Conclusion 

In this paper, we have considered a multicast channel, where a common data is transmitted from a 
source to several users. It is assumed that a minimum service must be provided for all the users. For this 
setup, we have optimized the average service received by a typical user in the network. Two scenarios are 
considered for the coverage constraint. In the case of hard coverage constraint, the minimum multicast 
requirement is stated in terms of the minimum rate (multicast outage capacity) received by all the users 
in a single transmission block. For small enough outage probabilities, it is shown that the capacity region 
is achieved by providing the required multicast rate in a single layer code, and designing an infinite-layer 
code as in [5], on top of it. In the case of soft coverage constraint, the multicast requirement is expressed 
in terms of the expected minimum rate received by all the users. An infinite layer superposition coding is 
shown to achieve the capacity region in this scenario. We have also proposed a suboptimal coding scheme 
for the MISO multicast channel. This scheme is shown to be asymptotically optimal, when the number 
of transmit antennas grows at least logarithmically with the number of users. 
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